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Abstract. Recently, strong, globally-ordered magnetic fields have been de- 
tected in some Herbig Ae and Be (HAeBe) stars, suggesting a possible evolu- 
tionary connection to main sequence magnetic chemically peculiar Ap and Bp 
stars. We have undertaken a detailed study of the binary system HD 72106, 
which contains a B9 magnetic primary and a HAeBe secondary, using the ES- 
PaDOnS spectropolarimeter mounted on the CFHT. A careful analysis of the 
very young primary reveals that it has an approximately dipolar magnetic field 
geometry, strong chemical peculiarities, and strong surface chemical abundance 
inhomogeneities. Thus the primary is very similar to an Ap/Bp star despite 
having completed less then 1.5% of its main sequence life, and possibly still 
being on the pre-main sequence. In contrast, a similar analysis of the secondary 
reveals solar chemical abundances and no magnetic field. 
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1. Introduction 

Strong, globally-ordered magnetic fields have been known to exist in Ap/Bp 
stars for quite some time. However, the origin and evolution of these magnetic 
fields, as well as their impact on the observed peculiar chemical abundances, 
remains unclear. Recently, magnetic fields have been discovered in some Herbig 
Ae and Be (HAeBe) stars (Donati et al., 1997; Hubrig et al., 2004; Wade et al., 
2005; Wade et al., 2007). HAeBe stars are pre-main sequence intermediate-mass 
stars and are the direct evolutionary progenitors of A and B stars. 
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The discovery of magnetic fields in HAeBe stars is of particular importance 
because it suggests a possible evolutionary link between Ap/Bp stars and mag- 
netic HAeBe stars. The observed magnetic fields in HAeBe stars appear to have 
similar strengths and geometries to those observed in Ap/Bp stars (Alecian et 
al., P27, these proceedings). Furthermore, magnetic fields are detected in HAeBe 
stars with approximately the same frequency as in main sequence A and B stars, 
5-10% (Wade et al., 2007). 

These recent discoveries prompted us to perform a detailed investigation of 
one particularly interesting case: the very young binary system HD 72106. 



2. HD 72106 

HD 72106 is a double star system with a 0.8" separation between the compo- 
nents. The brighter component (HD 72106A) was observed to possess a magnetic 
field by Wade et al. (2005), and the secondary was identified as being a HAeBe 
star by Vieira et al. (2003). 

The HD 72106 system was observed by Hipparcos, and solved as a binary, 
producing a "good quality" solution. Very large proper motions in both RA: 
—5.18 ± 1.08 mas yr -1 , and in DEC: 9.73 ± 1.36 mas yr _1 were observed, with 
no evidence for relative motion between the components (ESA, 1997). The Hip- 
parcos parallax for the system places the two stars at a distance of 288tg4 2 P c - 
The radial velocities we determine from our observations of both components 
are identical at 22 ± 1 kms -1 . Thus the stars are at approximately the same 
point in three dimensional space, and are moving together in three dimensions. 
This strongly suggesting that the system is in fact a true binary. 

Based on both Balmer line and metallic line fitting, we find T c g = 11000 ± 
1000 K and log g = 4.0 ± 0.5 for the primary, and T cff = 8750 ± 500 K and 
log g = 4.0 ± 0.5 for the secondary. 

The two stars can be placed on the H-R diagram, as seen in Fig. 1. When 
compared with the evolutionary tracks of Palla & Stahler (1993), we find that 
both stars sit near the zero-age main sequence (ZAMS). Based on the position of 
the secondary (which we infer to be a pre-main sequence star based on emission 
at Ha) we can constrain the age of the system to be between 3 and 13 Myr. At 
the younger end of this range, the primary would still have ~ | of its pre-main 
sequence lifetime remaining. However, at the older limit, the primary would 
already be on the main sequence, with a fractional main sequence age (r) of 
0.015. Thus the evolutionary status of the primary is unclear, although the star 
is certainly very young. 

In this investigation we used a set of 20 high-resolution Stokes I and V 
spectra obtained with ESPaDOnS, the high- resolution echelle spectropolarime- 
ter mounted on the CFHT. 
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Figure 1. The pre-main sequence H-R diagram of the HD 72106 system. Evolutionary 
track (solid lines) and isochrones (dashed lines), as well as the ZAMS and the birthline 
(for an accretion rate of 16~ 5 Mq yr" 1 ) from Palla & Stahler (1993) are shown. 



3. Magnetic Fields 

We applied the Least Squares Deconvolution (LSD) technique (Donati et a!., 
1997) to our observations of both HD 72106A & B. We find clear evidence for a 
magnetic field in nearly all our Stokes V profiles of HD 72106A, with maximum 
longitudinal field of 38 5 ±45 G, but we see no evidence for a magnetic field in any 
of our observations of HD 72106B (maximum longitudinal field of —51 ±55 G). 

In order to determine the magnetic field geometry of HD 72106, we need to 
know the rotation period of the star. We searched for periodicity in our lon- 
gitudinal magnetic field observations using a modified Lomb-Scargle method. 
We also applied this technique to both our Stokes / and V LSD profiles them- 
selves; searching each pixel in the line profile for periodic variability, and then 
averaging the results over all pixels. Additionally, the phasing of LSD profiles 
produced by candidate periods was investigated by eye, looking for variability 
physically consistent with rotation. The one period meeting all these criteria, 
0.63995 ± 0.00014 days, was adopted as the rotation period of the star. 

Using this rotation period, together with the inferred stellar radius and 
vsini, we determined the inclination of the rotation axis to the line-of-sight 
of HD 72106A to be i = 23 ± 11°. We modeled the magnetic field of HD 72106A 
by fitting the observed Stokes V LSD profiles with synthetic profiles computed 
by ZEEMAN2 (Landstreet, 1988; Wade et al., 2001). The mean Stokes I LSD 
profile was fit, then the phased set of Stokes V LSD profiles was fit using a 
dipole field model. From this procedure we derive the dipole polar strength 
Bp = 1300 ± 100 G and obliquity (3 — 60 ± 5°, with 3<r error bars. An analysis 
of the longitudinal field measurements shows that they are also consistent with 
this model. 
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Figure 2. Longitudinal magnetic field measurements, phased according to the rotation 
period of 0.63995 ± 0.00014 days, together with the best fit sinusoid for the data. 



4. Abundance Analysis 

The spectra of both components of HD 72106 were fit, using synthetic spectra 
calculated with the ZEEMAN2 spectrum synthesis code. Spectra were fit with 
a Levenberg-Marquardt \ 2 minimization method, and the results were exten- 
sively verified by eye. To approximate our magnetic observations, a 1 kG dipole 
magnetic field was included in the model of HD 72106A, and no field was used 
for HD 72106B. The final best fit values, reported in Fig. 3, are averages over 
seven spectral windows between 100 and 200 A in length, and the uncertain- 
ties are the standard deviations. When an abundance was based on only a few 
windows, a conservative estimate of the uncertainty was made by eye. 

For HD 72106A we find remarkably strong over-abundances of Cr, Fe and 
Nd. Additionally, Si is overabundant by ~ 0.8 dex while He appears to be 
underabundant by ~ 1.5 dex. We find usini = 41.0 ± 0.3 kms" 1 for this star. 
This pattern of strong peculiarities is distinctive of Bp stars. 

In contrast, for HD 72106B we find abundances that are consistent with 
solar, at least at the 2a level. Small departures from solar abundance for C and 
Sc may exist, however they are consistent with the normal scatter for an A star. 
We find v sin i = 53.9 ± 1.0 kms -1 and a microturbulcncc of 2.3 ± 0.6 kms . 

5. Doppler Imaging 

We performed Doppler Imaging of HD 72106A using the INVERS12 code (Kochukhov 
et al., 2004). LSD profiles were used for the Doppler Imaging, rather then in- 
dividual lines, due to the S/N of our spectra. LSD profiles for specific elements 
were constructed using line masks containing only lines of one element. 
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Figure 3. Mean best fit abundances, relative to solar, for both HD 72106A & B. Points 
marked with only an arrow indicate the value is an upper limit only. Solar abundances 
are from Grevesse et al.(2005). 



Surface abundance maps were calculated for Si, Ti, Cr, and Fe, using a 30° 
inclination angle, and are shown in Fig. 4. Strongly inhomogeneous distributions 
can be seen for all four elements. Ti, Cr, and Si have similar distributions. The 
distribution of Ti appears to differ slightly from that of Cr and Fe in the southern 
hemisphere, but the inversion is poorly constrained in that region. The large spot 
seen in the northern hemisphere near phase in all four images is close to the 
location of the positive magnetic pole. 



6. Conclusions 

HD 72106A is in many respects very much like an Ap/Bp star. It possess a 
predominantly dipolar magnetic field, strong chemical peculiarities, and strong 
surface abundance inhomogeneities. Despite this, the star is very young, having 
completed at most 1.5 % of its main sequence life and possibly still being on 
the pre-main sequence. In contrast, HD 72106B is chemically normal, and has 
no detectable magnetic field. 

The HD 72106 system provides evidence that chemical peculiarities can form 
at the ZAMS, and possibly even on the pre-main sequence. Additionally, it 
provides support for an evolutionary link between magnetic HAeBe stars and 
Ap/Bp stars. Finally this system raises questions about when in a star's evolu- 
tion chemical peculiar first appear, in particular whether peculiarities can arise 
during the pre-main sequence or only early in the main sequence. 



6 



CP. Folsom, G.A. Wade, et al. 



ip-0.0 <p-0.2 ip-OA <p-0.6 f-0.S 




Figure 4. Surface abundance maps of HD 72106A for Si, Ti, Cr and Fe. 
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